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Abstract. Extensions of the standard model of particle physics, in particular those based on string theory, often predict a 
new f/ ( 1 ) gauge symmetry in a hidden sector. The corresponding gauge boson, called hidden photon, naturally interacts with 
the ordinary photon via gauge kinetic mixing, leading to photon - hidden photon oscillations. In this framework, one expects 
photon disappearance as a function of the mass of the hidden photon and the mixing angle, loosely constrained from theory. 
Several experiments have been carried out or are planned to constrain the mass-mixing plane. 

In this contribution we derive new constraints on the hidden photon parameters, using very high energy y-rays detected 
from the Crab Nebula, whose broad-band spectral characteristics are well understood. The very high energy y-ray observations 
offer the possibility to provide bounds in a broad mass range at a previously unexplored energy and distance scale. Using 
existing data that were taken with several Cherenkov telescopes, we discuss our results in the context of current constraints 
and consider the possibilities of using astrophysical data to search for hidden photon signatures. 
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INTRODUCTION 

Typical extensions of the current standard model of par- 
ticle physics often contain extra U{1) gauge degrees of 
freedom. Many models based on string compactifications 
show that standard model particles are uncharged under 
the additional U{1) symmetries. Therefore these degrees 
of freedom belong to a "hidden sector", i.e., an experi- 
mentally so far unobserved set of fields uncharged' un- 
der the standard model gauge group. The existence of 
high mass particles charged under both the visible and 
hidden sectors (mediators) or gravity can produce nev- 
ertheless small interactions between the two sectors. As- 
suming that hidden sector particles are light^, note that 
current accelerator based experiments can be largely in- 
sensitive to the subtle effects of their existence. Therefore 
in the recent past, many high-precision experiments have 
been dedicated to search for new light particles (e.g. [2] 
or [3] for an overview) and many ideas for new experi- 
ments are under consideration [4] . 

Natural models contain at least one hidden photon [5] 
(or sometimes called paraphoton), i.e., the coiTespond- 
ing gauge boson to the new gauge group. Here we 



' In case of direct renormalizable couplings of the corresponding 
gauge boson to standard model matter precise measurements of the 
electroweak theory have shown that their masses must exceed a few 
hundred GeV [J ]. 

^ The masses typically considered belong to the sub-eV range. 



consider a minimal theory with just one^ U{l)k gauge 
group ([6], [5], [2], [7]), in addition to the normal elec- 
tromagnetic gauge. The most general low energy La- 
grangian allowed by the symmetries is 

+ 2^^2g,g^ (1) 

where F^v is the field strength tensor for the ordinary 
photon gauge field A^, defined by F^v = d^Ay — dyA^, 
and Bj^y the same tensor for the hidden photon field 
B^. The third term, also allowed by gauge invariance, 
coiTesponds to a non-diagonal kinetic term, the so-called 
kinetic mixing, where x is the mixing angle between 
photons and hidden photons. We assume x to be small. 
The last term describes a possible mass of the hidden 
photon. It can arise either from Higgs or Stiickelberg 
mechanisms, but in the former case the model suffers 
from additional constraints [8]. 

The non-zero kinetic mixing states that the and 
Bf^ fields are non-orthogonal. The transformation to an 
orthogonal basis, with canonical kinetic term, is done by 
the redefinition 

B^^S^-sinxA^ (2) 



' Of course theories with more than one additional U{1) gauge sym- 
metry exist, but they would be more cumbersome to handle [6]. 



Moreover, one observes that (1) now contains a 
non-diagonal mass term that mixes photons with 
hidden photons (find the relevant part of the redefined 
^ below, expanded 



1 



(5^5^ - IXS^A^, + X^A^A^,) (3) 



leading to vacuum y-y^ oscillations (if hidden photons are 
not massless). Here is the quantum of the field 5^, that 
being orthogonal to the photon is completely sterile with 
respect to electromagnetic interactions. Note also that the 
diagonalization causes a multiplicative renormalization 
of the electric charge. 

It is interesting to note that the oscillation effect is 
completely analogous to the phenomenology of neutrino 
oscillations. Thus the oscillations of photons open the 
possibility to search for hidden photons via, e.g., "light 
shining through a wall" (LSW) experiments ([6], [1], 
[2]). 

In order to compute the oscillation probability at a 
distance L, one has to solve the equations of motion for 
the Lagragian to find the propagation eigenstates (see e.g. 
[2] for a short review). One obtains the following result 
for the oscillation probability (in natural units): 



(4) 



where E stands for the energy. Hence the oscillation 
length is observed to 



^osc — 



ju2 l^TeV ) VlO-^eV 



kpc. (5) 



The values of the mixing parameters of hidden pho- 
tons % and ju are unspecified from theory and experimen- 
tally unknown because these particles have not been de- 
tected until today. However, there already exists a broad 
range of restrictions on these parameters provided by 
several experiments. The stronger constraints originate 
from tests of the Coulomb law, CMB measurements, 
LSW experiments, and searches of y^'s radiated from the 
Sun (see e.g. [8] and references therein). Find a compos- 
ite plot of the current bounds in [8], Fig. 1 . 

Considering the experiments, it turns out that the 
bounds were obtained using energies up to ^(lOOGeV) 
(LEP) and distances up to ^(1 AU) (solar searches), so 
it's worth mentioning that no constraints exist using the 
very high energy range (i^(TeV)) and distances with 
^(kpc). In this work we will provide such (astronomi- 
cal) bounds using very high energy (VHE, E > lOOGeV) 
y-rays from the Crab Nebula. 

The Crab Nebula (Messier 1) is a remnant of a super- 
nova explosion that occurred in the year AD 1054. A pul- 
sar (PSR B0531H-21) is located near the geometrical cen- 
ter of the nebula. Today, the remnant is observed to be of 



the plerionic type at a distance of dc = (1 .93 ± 0.1 l)kpc 
from Earth 19], with a bright continuum emission from 
radio to very high energy y-rays, peaked in the near- 
infrared and optical range. The whole emission is pre- 
dominantly produced by non-thermal processes, mainly 
by synchrotron and inverse Compton interactions of en- 
ergetic electrons. The mechanism producing the VHE 
spectrum is inverse Compton scattering of accelerated 
electrons (up to PeV energies) on different low energy 
seed photon fields, dominated by the synchroton field 
(see e.g. Figure 10 in [10]). 

The broad-band VHE spectrum can be parametrised 
by (see [10] for the coefficients pi) 

where v/v = is the differential energy-flux of the 
Crab nebula. 

Upper Umits measured for the diameter of the nebula 
in file VHE regime are Cfc < 2' at £ < lOTeV and Ofc < 3' 
at £ > lOTeV [10]. However, considering models of the 
VHE emission of the nebula (e.g. [10]) leads to diameters 
of ^(arcsec), which we assume here. 

Before we explain the method of giving new con- 
straints, there remains one subtlety we have to men- 
tion briefly. The oscillation probability (4) is calculated 
under the assumption that the photons can be repre- 
sented as plane waves. Considering production and de- 
tection processes, this assumption does not hold under 
normal circumstances. Rather, (4) must be calculated us- 
ing quantum mechanical wave packets having a coher- 
ence length D. When working with wave packets of dif- 
ferent masses mi ^ m2, the velocities would differ by 
a factor Aj3 ~ Am^/E-^, Ani^ = m\— m^, so they sepa- 
rate by LA/3 after traveling a distance L. Thus osciUations 
freeze out, if 

D 
A/3 

(see [11]). In our case, Arn^ = ji^. Detailed quantum 
mechanical calculations with wave packets reveal (see 
[12]), that for relativistic particles 
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where Cx = y o^p + o^^, with cr^p, Oxd the spatial un- 
certainties of the production and the detection process, 
respectively. To get oscillations, three conditions have to 
be satisfied: (a) Lcoh > <ic ~ 2kpc, (b) Lose < Lcoh, and 
(c) Ox «. Lose- 



For the production process one finds Oxp calculating 
the interaction length for electrons in the nebula, which 
interact via inverse Compton scattering with various seed 
photons. A lower limit is found assuming the Thomson 
regime, ye <^ mi.c^, with 7 the Lorentz-factor of the 
electron, e the energy of the seed photon, nie the electron 
rest mass, and c the speed of light in vacuum. In this case 
the interaction length is given by 

= c7T^de'«b(e'), (9) 

where at is the Thomson cross section and nb(e) = g| is 
the differential (seed-) photon density. For a rough lower 
limit, approximating the main seed field (synchrotron) 
in [10] and using the extension of the nebula it reveals 
A > 1 kpc. 

Comparing X to the radius of the nebula in the VHE 
regime rc <C 1 pc, it should be clear that rc is the right 
estimation for Gxp, because ^ A. 

The spatial uncertainty of the detection process, Gxd, 
does not contribute. Considering the detection process in 
detail one finds that a^o ^ <7rP (ZechUn, et al., in prep.). 

Taking these results into account, one observes that the 
conditions given above hold for masses jU < 10-5 eV. 

METHOD 

As mentioned above, this work investigates the very high 
energy data available from the Crab nebula. The VHE 
regime is mainly covered by ground-based detection 
techniques, due to the fact that non-thermal 7-ray sources 
typically produce power-law type spectra in which the 
flux drops with increasing energy. As a consequence 
large effective detection areas of ^(10^ m^) are required, 
sufficient to compensate for small photon fluxes. The 
data used here were taken with (stereoscopic) imaging air 
Cherenkov telescopes (lACTs), especially by HEGRA 
[10], H.E.S.S. [13], MAGIC [14], and the Whipple 10m 
telescope [15]. 

If hidden photons of mass ji exist, the energy de- 
pendent oscillation probability (4) influences the observ- 
able spectrum. Some photons will convert to hidden pho- 
tons during their propagation to Earth which are not de- 
tectable because they initialize no air-shower Due to 
the broad-band VHE spectrum that shows no intrinsic 
absorption lines in this region these spectral signatures 
should be measurable or could be used to constrain x 
and jj.. The data and the expected spectral signature are 
shown in Fig. 1. 

The differential spectra are measured as functions of 
energy-bins (see e.g. [10]), where E denotes the geomet- 
ric mean energy of a bin with width AE. Thus, one has to 
average (4) (at the Crab distance dc) over the bin-size of 
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FIGURE 1. The spectral data of four experiments is shown 
(references in the text). In addition, the signature of hidden 
photons = (v/v) • (1 — Py^y^{dc)) for specific parameters 
is demonstrated in comparison with the inverse Compton flux 
v/v (6). Note that the data is not scaled (for further details 
consider the following section). 



the energy-bin centered on E, giving 

Py^yXEAE) = ^ J dEPy^yXE.d,) (10) 

A£(£) 

Therefore, the predicted energy-flux observed is 

y (F, A£) = ( v/v ) I • ( 1 - Py^y, (E,AE)) (11) 

for the parameters x and jj., where the inverse Compton 
flux v/v from (6) was used. 

Hence one has to fit the spectral signature (1 1) to the 
data. This can be done applying a goodness-of-fit test, 
here the method of least squares is used (see e.g. [16]). 
In our case, Xisq is given by the expression 

where the sum runs over all data points, given by 
(£■,■,}',■, a,), yi stands for the measured energy-flux at 
energy £, with a statistical error a,-. The method of 
least squares can be applied considering arbitrary con- 
fidence levels. The following results are calculated us- 
ing the 68.3% confidence level which is constrained by 
^X\sq = 1 (for one fit parameter). 

Solving this numerically under the conditions ex- 
plained above for every allowed mass ji (< 10-5 eV, see 
above) one gets a fit value Xtit for X for every jj. (Xat de- 
pends on the desired confidence level). Interpreting this 
in the disappearance approach it is clear that one can ex- 
clude all values x > Zfit- 



RESULTS 

Clearly, comparing the data of the different telescopes 
(see Fig. 1) it is worth mentioning that the results of 
the experiments differ among each other. Physically and 
within the experimental errors, all experiments must 
measure the same flux at a constant energy E. We choose 
to renormalize the energy scale of the instruments within 
their respective systematic uncertainties to avoid smear- 
ing of signatures. The scaling factor can be found by a fit 
of (6) on every data set. 

TABLE 1. To avoid smearing of signatures the energy axis 
of every experiment has been rescaled by E' = .v ■ E. Values 
X ^ Xfa can be excluded. To check the goodness of the fit, the 
minimum reduced chi-squared is given for fi = 10~* eV. 



Telescope 


Scalefactor s 


X&t 






HEGRA 
H.E.S.S. 
MAGIC 
Whipple 


1.000 
0.921 
0.968 
0.944 


0.0346 
0.0800 
0.1165 
0.1036 


1.43(15) 
4.05(7) 
0.82(9) 
1.21(7) 


combined 




0.0343 


1.60(41) 



For energies above lOTeV, the re-scaled spectra of 
H.E.S.S. and HEGRA still deviate from each other. Since 
the H.E.S.S. data are not consistent with the model con- 
sidered (see (6)), we choose to ignore the data for ener- 
gies above lOTeV. 

Applying the method described above and using the 
HEGRA, H.E.S.S., MAGIC, and Whipple data we can 
conclude with a value Xm — 0.0343 (the result converges 
to this value with increasing mass). The result is shown in 
Fig. 2. Results for separate experiments can be found in 
Tab. 1 . To test the goodness of the fit, the minimum of the 




0.01 



H[eV] 



reduced chi-squared value, defined by Zred(") /tisq/"' 
is given in the table for a specific mass fi, where n is the 
number of degrees of freedom. 

Due to experimental errors on the distance dc and the 
energy £, of each data point, the error on the mass ji of 
every point excluded from the mass-mixing plane can be 
approximated to be 9 %. 



CONCLUSIONS 

Comparing our results to the constraints given in [8], 
Fig. 1, the limits obtained from measuring deviations 
from the Coulomb law are better. But the bounds given 
here are the best constraints on the hidden photon pa- 
rameters using oscillation effects of photons directly. For 
the first time, we got astronomical Umits considering new 
energy and distance ranges. 
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